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ABSTRACT 

We present the discovery by the WASP-South survey of three planets transiting moderately 
bright stars (V ~ 11). WASP-120b is a massive (d.ShAfjup) planet in a 3.6-day orbit that we 
find likely to be eccentric (e = O.OSQtQ g^g) around an F5 star. WASP-122 b is a hot-Jupiter 
(1.28Mjup, 1.74i?jup) in a 1.7-day orbit about a G4 star. Our predicted transit depth variation 
caused by the atmosphere of WASP-122 b suggests it is well suited to characterisation. WASP- 
123 b is a hot-Jupiter (0.90Mjup, 1.32i?jup) in a 3.0-day orbit around an old 7 Gyr) G5 
star. 

Subject headings: Planetary systems — stars: individual (WASP-120,WASP-122WASP-123) 


1. Introduction 

The Wide Angle Search for Planets (WASP) 
survey is a prolific contributor to the field of ex¬ 
oplanet science having published the discovery of 
104 planets to date. Our effective magnitude range 
of 9 < V < 13 means that WASP systems are con- 
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ducive to further study. Examples from the ex¬ 
tremes of this range are the bright WASP-33 (V 
= 8.3; ICollier C anieron et al. ioifil l and WASP- 
18 (V = 9.3: lHellier et aLl2009) and the relativel y 
dim WASP-112 ( V = 13.3; Anderson et al.ll20l3l . 


Here we present the discovery of: WASP-120 b, 
a system with a star showing variable activity 
and a possibly eccentric planet orbit, WASP-122 b, 
which offers a good opportunity for atmospheric 
study, and WASP-123 b, which orbits an old star, 
^ 7 Gyr. 


2. Observations 


The transits of these planets were discovered 
in photometry gathered from the WASP-South 
installation hosted by the South African Astro¬ 
nomical Observatory. The WASP-South instru¬ 
ment is an array of 8 cameras using 200mm 
//1.8 lenses to survey the sky at a cadence of 
^ 10 minutes. For more information on the 


WASP instrument, see lPollacco et al.l (|2006ll . The 
data were pr ocessed and searched for tran sits as 
described in Gollier Cameron et al.l ( 2006 ) with 
candidate selection fo l lowing the procedure in 
Collier Cameron et al.l (l20n7al ). Details of obser¬ 
vations for each star in this paper can be found 
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Fig. 1.— Discovery data for WASP-120 b. 

Top panel: Phase folded WASP photometry for 
WASP-120. Middle panel: WASP discovery pho¬ 
tometry (grey), TRAPPIST (blue) and Euler- 
Cam (green) follow-up photometry with our tran¬ 
sit model over-plotted. The meridian flip in 
the TRAPPIST data has been corrected for and 
marked with a vertical dashed line. All photo¬ 
metric data have been binned with a duration of 
2 minutes for clarity. Bottom panel: CORALIE 2 
radial velocity data from before (red circles) and 
after (blue triangles) the upgrade, over-plotted 
with the circular (black-solid) and eccentric (grey- 
dashed) solutions. 


Fig. 2.— Discovery data for WASP-122 b. 
Top panel: Phase folded WASP photometry for 
WASP-122. Middle panel: WASP discovery pho¬ 
tometry (grey), TRAPPIST (blue) and Euler- 
Cam (green) follow-up photometry with our tran¬ 
sit model over-plotted. The meridian flip in 
the TRAPPIST data has been corrected for and 
marked with a vertical dashed line. All photomet¬ 
ric data have been binned with a duration of 2 
minutes for clarity. Bottom panel: CORALIE ra¬ 
dial velocity data, over-plotted with our circular 
solution. 
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Fig. 3.— Discovery data for WASP-123 b. Cap¬ 
tion as for Fig. [2j 
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Fig. 4.— Bisector spans plotted against relative 
radial velocities for WASP-120 (top), WASP-122 
(middle) and WASP-123 (bottom) showing no cor¬ 
relation. Solid lines are results of least-squares fits 
to the data, dashed lines are the Icr uncertainties 
of the fits. Date of observation is denoted by point 
colour. The increased dispersion of points seen in 
more recent data for WASP-120 is attributed to 
an increase in stellar activity. 
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in Table [TJ The phase-folded WASP data are dis¬ 
played in the top panels of Figs. ,5UI and m We 
used the method of iMaxted et alJ (j2nilh to search 
the WASP photometry for modulations caused by 
star spots. We detected no rotational modulation 
above 2mmag which suggests that the hosts are 
inactive. 

We obtained spectra of the three stars with the 
CORALIE spectrograph on the 1.2-m Swiss tele¬ 
scope as outlined in Tabled] We used these data to 
measure radial velocity (RV) variations and con¬ 
firm the planetary nature of the candidates (Table 
[21 bottom panel of Figs. HI [2] and [3]). We ob¬ 
tained 9 of the WASP-120 spectra after the spec¬ 
trograph was upgraded in November 2014. The 
lack of correlation between the bisector spans and 
RVs (Fig. |4|) indicate that the RV variations are 
not a r esult of blend e d ecli psing binaries. For ex¬ 
ample, ISantos et al. ( 2002 1 found a brown dwarf 
mass companion that produces a correlation be¬ 
tween the RVs and bisector spans with a gradient 
of 0.67. The largest gradient from amongst our 
planets is an order of magnitude smaller and not 
significant; 0.06 ± 0.07 for WASP-120. While we 
cannot strictly rule out the case of blended planet 
hosting stars we can eliminate brown dwarf blends 
and more massive objects. 

We acquired the follow-up photometry needed 
to accurately determine the system parameters 
from the 0.6-m TRAPP IST telescope dQllon et aF 


201ll) and FulerCam ( Lendl et al. 2012 1 on the 
Swiss telescope at La Silla, Chile. The TRAP¬ 
PIST telescope’s equatorial mount requires a 
meridian flip when the target culminates dur¬ 
ing an observation. These occurred at BJD = 
2456609.725 during the transit of WASP-120 on 
2013 Nov 12 and at BJD = 2456644.758 during 
the transit of WASP-122 on 2013 Dec 17. We ac¬ 
count for any offsets introduced by treating them 
as two separate datasets during our analysis. The 
photometric data are presented in Table |3| This 
follow up revealed a star within 2.2” of WASP-120 
which is 4.35 ± 0.02 magnitudes fainter in the 1 
band and 3.89 ± 0.02 magnitudes fainter in the z 
band. CORALlE’s fibres are 2” in diameter and 
the RVs were obtained in good seeing, so the star 
is sufficiently distant that it did not contaminate 
the observations and thus could not cause a false 
positive. 


Table 1; Observations of WASP-120, WASP-122 
and WASP-123 


Date 

Source 

N.Obs. / 
Filter 

Comment 

WASP-120 

2006 Aug-2012 Jan 

WASP-South 

27 079 


2013 Sep-2015 Mar 

CORALIE 

29 


2013 Nov 12 

TRAPPIST 

I+Z 

Meridian flip 

2014 Aug 17 

EulerCam 

Ic 


2014 Dec 07 

EulerCam 

Ic 


WASP-122 

2011 Oct-2012 Mar 

WASP-South 

4 834 


2013 Nov-2014 Oct 

CORALIE 

17 


2013 Dec 17 

TRAPPIST 

z 

Meridian flip 

2014 Jan 15 

TRAPPIST 

z 

Slight cloud 

2014 Jan 15 

EulerCam 

r 

Slight cloud 

2014 Jan 22 

EulerCam 

r 

Slight cloud 

WASP-123 

2006 May-2012 Jun 

WASP-South 

13 267 


2013 Sep-2014 Aug 

CORALIE 

20 


2014 Sep 06 

EulerCam 

Ic 


2014 Sep 06 

TRAPPIST 

z 


2014 Sep 09 

TRAPPIST 

z 



Table 2: Radial velocity data from CORALIE. 


HJD 

- 2 540 000 

RV 

(km s“^) 

Error 
(km s“^) 

BS 

(km s“^) 

Target 

Name 

6552.902673 

19.30305 

0.05043 

-0.27331 

WASP-120 

6572.735422 

20.41666 

0.03569 

0.03952 

WASP-120 

6573.843533 

19.62205 

0.03173 

-0.01785 

WASP-120 

6871.746191 

16.99747 

0.00673 

-0.01808 

WASP-123 


Note. —Data available in this format at ADS. The data are provided 
to the full precision used in our calculations but times are only accurate 
to a few seconds at best. 
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Table 3: Follow-up photometry from TRAPPIST and EulerCam. 


HJDutc 
- 2 540 000 

Norm. 

Flux 

Error 

AX 

Position 

AY 

Position 

Airmass 

Target 

FWHM 

Sky Bkg. 
(Counts) 

Exp. 
Time (s) 

Target Name 

Instrument 

Band 

6887.720753 

0.990347 

0.002159 

-1.91 

-2.21 

2.64 

11.12 

71.30 

50.00 

WASP-120 

EulerCam 

IC 

6887.721534 

0.996407 

0.002145 

0.25 

-1.40 

2.62 

12.18 

70.00 

50.00 

WASP-120 

EulerCam 

IC 

6887.722324 

0.994713 

0.002114 

0.06 

-1.04 

2.60 

10.85 

68.82 

50.00 

WASP-120 

EulerCam 

IC 

6910.766640 

0.992689 

0.008096 

-0.92 

-0.64 

2.91 

4.08 

293.86 

13.00 

WASP-123 

TRAPPIST 

z 


Note. —Data available in this format at ADS. The data are provided to the full precision used in our calculations but times 
are only accurate to a few seconds at best. 


3. Analysis 


3.1. Stellar Parameters 

We determined the atmospheric parameters 
of each host star by analysing the co-added 
CORALIE spectra after correcting them for shifts 
due to the radial motion of the star using the mea¬ 
sured RVs. O ur spectral a n alysis followed proce¬ 


dures given in iDovle et al.l (12013^ . For each star 


we obtained the effective temperature, Tes, using 
the Ha line, log 5 from the Na D and Mg b lines 
and iron abundances from the analysis of equiva¬ 
lent width measurements of several unblended Fe i 
lines. We found the projected rotation velocity, 
V sini, by fitting the profiles of the Fe i lines after 
convolving with the instrumental resolution (H 
= 55 000) and a macro turbulent ve l ocity adopted 
from the calibration of Dovle et al. ( 2014ll . 


3.2. System Parameters 

We used a Markov chain Monte Carlo (MCMC) 
code to determine the system parameters using the 
discovery an d follow-up photometry with RV s as 


described by Collier Cameron et al. ( 2007b[l and 


Anderson et al.l (1201511 . 

For each system we modelled our transit 


l ightcu rves using the formulation of iMandel &: Ago! 
(|2002ll and accounted for limb-da rkening using 
the fo ur-parameter non-linear law of Claret ( 2000l 
200411 . The photometric bands and limb-darkening 
coefficients used in the lightcurve models are de¬ 
tailed in Table m 

We used bagemass ( Maxted et al.l I 2 OI 5 II to 
compare ps , determined from the transit lightcurves 


coupled with the spectroscopic values of [Fe/H] 
and Teff, to stellar models in order to estimate the 
mass of the star. BAGEMASS also gives an estimate 
of the age of the system. 


To calculate the distan ce we use the a p paren t 
Kg band magnitude from ISkrutskie etT] (1200611 . 
the radius of the star from Tableland the angular 
diameter of the star based on the calibration of 
K-band surfa ce brightness - e ffectiv e temperature 
relation from Kervella et al.l ( 2004h . We assume 
that interstellar reddening is negligible and that 
K= Kg -f 0.044. 


The free parameters in our MCMC analysis 
were Tq, P, {Rp/RsY, Tu, b, Ki, 7 , [Fe/H] and 
Tld- Here Tg is the epoch of mid-transit, P, is 
the orbital period, {Rp/RsY is the planet-to-star 
area ratio, T 14 is the total transit duration, b is 
the impact parameter of the planet’s path across 
the stellar disc, Ki is the reflex velocity semi¬ 
amplitude, 7 is the systemic velocity, [Fe/H] is the 
stellar metallicity and Tld is the limb-darkening 
temperature. Tld and [Fe/H] were constrained 
by the spectroscopic values of Tes and [Fe/H]. 
Tld was used by the MCMC to interpolate limb- 
darkening coefficients at each step from the Claret 
limb-darkening tables for the appropriate photo¬ 
metric band of each lightcurve. At each step of 
our MCMC these values were perturbed by a small 
random value and the of the model based on 
the new values was calculated. If this lead to 
a lower the step was accepted while a larger 
value would be accepted with a probability pro¬ 
portional to exp(—Ay^/2). Our final values were 
calculated from the medians of the posterior distri¬ 
butions with uncertainties corresponding to the Ict 
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Table 4: Limb-darkening parameters extrapolated using the Tld resulting from each analysis. 


Planet 

Instrument 

Instrument Band 

Claret band 

ai 

a2 

as 

a4 

WASP-120 

WASP 

Broadband (400-700 nm) 

Cousins R 

0.136 

1.286 

-1.188 

0.391 


TRAPPIST 

I+Z 

Sloan z 

0.221 

0.827 

-0.760 

0.226 


EulerCam 

Cousins I 

Cousins I 

0.200 

0.951 

-0.863 

0.263 

WASP-122 

WASP 

Broadband (400-700 nm) 

Cousins R 

0.717 

-0.503 

1.076 

-0.515 


TRAPPIST 

z 

Sloan z 

0.799 

-0.743 

1.095 

-0.492 


EulerCam 

Gunn R 

Cousins R 

0.717 

-0.503 

1.076 

-0.515 

WASP-123 

WASP 

Broadband (400-700 nm) 

Cousins R 

0.683 

-0.405 

0.957 

-0.473 


TRAPPIST 

z 

Sloan z 

0.766 

-0.664 

1.010 

-0.462 


EulerCam 

Cousins I 

Cousins I 

0.763 

-0.639 

1.059 

-0.491 


conhdence intervals. When we allow the MCMC 
to explore eccentric solutions we fit ^/ecosuj and 
^/e sin oj to ensure a uniform probability distribu¬ 
tion. Our results for each system are in the lower 
part of Table [5] and corner plots of the jump pa¬ 
rameter posterior distributions of each analysis in 
Figures [5l [S] and [T] 

We checked for trends in the derived tran¬ 
sit depths with respect to the colour of the ob¬ 
servational band for each star by running each 
lightcurve through our MCMC separately. The 
depths for WASP-122 and WASP-123 agree to 
within 1 cr of the depth derived from the combined 
analysis. The depths from this analysis of the two 
complete lightcurves of WASP-120 show a slight 
difference of (1.2 ± 0.4) x 10“^ which could be ac¬ 
counted for by a low level of inherent stellar vari¬ 
ability in either the host star or faint, nearby com¬ 
panion. 


4. WASP-120 


WASP-120 b is a 4.85-Mjup, 1.73-i?jup planet 
orbiting a moderately bright {V = 11.0) F5 star. 
The effective temper ature of WASP-120 p laces it 
in the lithium gap (|B5hm-Vitensel l2004ll . so we 
cannot estimate the age of this star based on 
the lithium abundance. Using the star’s Tycho 
B—V colour and rotation period from its V sini 
and radius from o ur MCMC we use gyrochronol- 
ogy calibration of BarnesI (2007) to estimate an 
age of 0.7 ± 0.6 Gyr. 


For compari s on, the cal¬ 


ibration of iMamaiek fc Hillenbrandl (|2008l) gives 
1.0 ± 1.8 Gyr. We ca nnot apply the calibration of 
Meibom et al. ( 20091) as the star’s colour results in 
a term requiring the logarithm of a negative value. 
Using BAGEMASS we hnd an age of 2.6 ± 0.5 Gyr 
which is consistent with that of the Mamajek & 


Hillenbrand calibration. 

The FWHM of the lines in the spectra and the 
bisector spans show more scatter in the later data, 
after the CORALIE upgrade (Fig. [8]), suggesting 
that the star may have become more active, and 
therefore have variable activity like the Sun. It is 
unlikely that the increased scatter is caused by the 
change to CORALIE since datasets on other stars 
don’t show an increased s catter (e.g. recent RV 
data taken 
20161 ) . 


of WASP-47; Neveu-VanMalle et al 


Due to the upgrade to CORALIE we parti¬ 
tioned the data into two sets. We added jitter of 
5.1±0.2 ms“^ to the older data and 6.3±0.8 ms“^ 
to the newer data. These values were adopted such 
that both datasets gave reduced values of one 
compared to a circular-orbit solution, and are in 
keeping with j itter determined for similar stars by 


WrightI ( 2 OO 5 I) . The uncertainties were estimated 
using a jackknife resampling. The difference be¬ 
tween the jitter values of 0.8 ± 0.8 is consistent 
with zero, thus the jitter values are consistent to 
Icr. 

Our resulting orbital solution had an eccentric¬ 
ity of 0.057to)oi8' This is signific antly non-zero at 
3.3g, while a Lucy-Sweeney test (jLncv fc Sweeney 


19711 ) gives a probability of only 0.1% that the or¬ 


bit is circular. We note, though, that this result 
is somewhat dependent on the jitter values used. 
Fitting with no jitter increased the eccentricity to 
0.068 ±0.014. This led to a higher apparent sig¬ 
nificance (4.9 (t) though the Lucy-Sweeney test no 
longer excluded the circular solution since neither 
solution is a good fit to the data. We thus adopt 
the value of eccentricity with jitter added but re¬ 
gard this as needing confirmation. 

More massive planets, such as WASP-120 b at 
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Fig. 5.— Plots showing the relationships of various jump parameters To, P, {Rp/Rs)^, Fi 4 , b, Ki, 7 , [Fe/H], 
Tld, y/ecosuj and y/esinuj from the analysis of WASP-120. All distributions have h ad the mean value from 
Table [S] subtracted. Plots were prepared with a modifed version of trangle.py by iForeman-Mackev et al.l 
(I2OI4II 
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Fig. 6.— Caption as for Fig. [S] As we enforced a circular orbit we did not fit ^cos lu and -yesino;. 
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Fig. 7.— Caption as for Fig. |6] 
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5.0-Mjup, often have eccentric orbits (Fig ITOl) . 
though it is unclear if the correlation is due to 
a real phenomenon or observation bias as the ec¬ 
centricities of more massive planets are easier to 
detect from RVs. A conclusive determination of 
this system’s eccentricity would come from obser¬ 
vation of an occultation which we expect to be 
delayed by 2.6 ± 1.1 hours if our eccentricity value 
is accurate. Based on the equilibrium tempera¬ 
ture of the planet, we estimate occultation depths 
in the Spitzer 3.6 fim and 4.5 fim bands of ap¬ 
proximately 780 and 960 ppm, respectively . Re- 


cent observations with Spitzer I Zhao et al. 20141: 


IPeming et al.l l2015h show detecting such an oc¬ 
cultation is easily achievable. 

The faint star close to WASP-120 may pro¬ 
vide evidence in support of the planet having un¬ 
dergone high-eccentricity mi CTation due to Koza i- 
Lindov cycles. According to IPotter et al.l ( 2008h . 
the companion’s colour is consistent with a 0.6 
Mq, K9 star at the same distance as WASP-120. 
Assuming this to be the case, the on sky separation 
gives a mini mum separation of 950 AU. U sing the 
equations of Fabrvckv fc Tremaine! ( 2007f ). a star 
of this mass could induce Kozai-Lindov cycles if 
the planet’s original orbital distance was 14.5 AU 
or greater. As the periodogram of the RVs for 
WASP-120 show no significant peaks beyond that 
of the planet (Fig. |9|) we checked for the presence 
of an additional, long period, object in the system 
by fitting a linear trend to the residuals of the 
RVs. The result was an RV drift, 7 , of (84 ± 73) 
ms“^ vr~ ^ which is cons i stent with zero at ^ 1.2tT. 
Following iMontet et al. ( 2014 ), for a planet on a 
circluar orbit: 


7 = (6.57ms ^ yr 


M2 


M 


Jup 


) (^) 


Therefore, if there is an additional object in the 
system it has a mass, M 2 , semi-major axis 02 , and 
inclination, 12 , such that, M 2 sin* 2/02 ^ 0.51Mjup 
AU-2. 

Notable examples of massive planets with con- 
fidently detected ecc entricities are HAT-P-16b 


( Buchhave et al.N201Cir). HAT-P-21b (|Bakos et al 


201 ih WASP^b (I.Toshi et al.ll2009l) and WASP- 
89 b ( Hellier et al. 2015 ). All of these are in sub 


7-day orbits with mas ses > 4 M.t 


iini_ 


table are HAT-P-20b ( Bakos et al.l 12011 ). which 


Also no- 


has the smallest eccentricity of the group and 
Kepler-14 b with the longest o rbital period at 


6.79-days (jBuchhave et al.l 120111) . Like WASP- 


120, 3 of these 6 systems are known to have other 
stars nearby; HAT-P - 16, WASP14 and Kepler-14 


( Wdllert et al. l2015t Wollert fc Brandnei 2015t 


Ngo et al. 2015[ Buchhave et al.l 201lh However 


this sample is too small to draw conclusions about 
a link between orbital eccentricity and the pres¬ 
ence of a stellar-mass neighbour. 

It is thought that stars with effective tempera¬ 
tures cooler than 6200 K have convective envelopes 
which enhance or bit circularisation/re-alignment 


Winn et al. ( 2010[) . This could erase any correla¬ 


tion between the type of orbit a planet is in and 
the presence of a further companion. 

In a study on th e prevalenc e of m ultiple stars 
in planetary systems iNgo et ah ( 2015[) found that, 
of their sample of hot host stars {TeS > 6200 AT) 
with evidence of misaligned or eccentric planet or¬ 
bits 59% ± 17% had companions while 83% ± 14% 
of their well-aligne d/circular o rbit sample had 
companions. When iNgo et al.l ( 2015 ) considered 
just the spin-orbit alignment of these systems, 
73% ± 15% of misaligned systems had compan¬ 
ions as opposed to 53% ± 14% of well aligned sys¬ 
tems. They concluded that there is no evidence 
for a link between multiplicity and orbital eccen¬ 
tricity/misalignment, though so far the sample is 
just 18 hot stars, fewer when just those systems 
with measured spin-orbit alignments are used. 

The typical time-scale of orbital circularisation 
is expected to be shorter than that of tidal re¬ 
alignment meaning that observations of spin-orbit 
misalignment may provide a better record of mi¬ 
gration pathway than eccentricity. Therefore, in¬ 
creasing the number of hot-host star systems with 
measured spin-orbit alignment that have been 
evaluated for stellar multiplicity could change the 
current picture. 

5. WASP-122 

WASP-122 b is a 1.28Mjup, 1.74i?jup planet 
orbiting a moderately bright (V = 11.0), metal 
rich ([Fe/H] = +0.32 ± 0.09), G4 star. WASP- 
122 is depleted in lithium (lo g A(Lf) < 1 . 0 ) and 


so m ust be several Gyr old (|Sestito fc Randich 
20051) . Using the star’s colour and rotational pe¬ 
riod from its Usinf, we calculate a gyrochrono- 
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Fig. 8.— Full width at half-maximum (FWHM) and bisector spans of spectra as a function of time for 
WASP-120. The increase in scatter in more recent data is taken as an indication the star may be entering a 
phase of increased activity. Red circles are data taken before the CORALIE upgrade and blue triangles are 
those taken after, the black dotted lines denote the date CORALIE was back on sky after the upgrade. The 
data have all had the mean of their distribution subtracted before plotting. 
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Fig. 9.— Periodogram for the RV data of the WASP-120 system. The only significant peak is that corre¬ 
sponding to the planetary period. 
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logical ages of 2.3 ± 1.4 Gyr ( 

BarnesI 2007). 2.8 ± 

1.4 (Mamaiek & Hillenbrandl 

20o 4 and 2.9 -b 4.8 

(Meibom et al. 2009|). Using BAGEMASS we find 


two possible solutions. Approximately 75% of the 
Markov chain output by BAGEMASS favour a mass 
of 1.24 ± 0.04 Mq and an age of 5.11 ± 0.80 Gyr. 
The other 25% of the output prefer a solution 
giving a mass of 1.10 ± 0.03 Mq and an age of 
8.67 ± 1.05 Gyr. The favoured, younger, higher- 
mass solution is a better match to, though still 
older than, the gyrochronological ages. 

WASP-122 b presents a good target for atmo¬ 
spheric characterisation via transmission spec¬ 
troscopy. Assuming the atmosphere is isothermal 
and adequately described as an ideal gas we can 
calculate the atmospheric scale height, H, using: 




( 2 ) 


Here, k is Boltzmann’s constant, Teq is the plan¬ 
etary equilibrium temperature, g is the planetary 
surface gravity and /i is the mean molecular mass 
of the atmosphere. We can use H to predict the 
transit depth variation due the addition of one 
atmospheric scale height to the planetary radius. 
In the case of WASP-122 b this is 142 ppm. The 
same calculation for the well-studie d HD209458b 
yields a variation of « 200 ppm. iDeming et al 


( 2013ll found evidence for water absorption on this 
scale in HD209458b. While WASP-122 is dim by 
comparison to HD209458 constraints have been 
put on the atmospheric compositions of planets 
with similarly bright hosts. For example, stud¬ 
ies of WASP-12 (V = 11.6, ) show ev idence of 
aerosols and a lack of TiO ( Sing et ^ 12013 1 as 
well as pla cing constraints on th e G/0 ratio of 
the planet (jKreidberg et al.ll2015r) which has been 
suggested may be an indicator of formation envi¬ 
ronment. We predict occultation depths in 3.6 gim 
and 4.5 fim Spitzer bands of 2100 and 2500 ppm 
respectively. In the K-band we predict a depth of 
1000 ppm. Similar K-band depths ha ve been de¬ 


tecte d, for example that of WASP-lOb (I Cruz et al 


2015ll . making ground-based follow up possible. 
Such observations of WASP-122 b stand to shed 
light on our understanding of atmospheric albedo 
and opacity sources as well as its formation his¬ 
tory. 


6. WASP-123 


WASP-123 b is a 0.90-Mjup, 1.32-i?jup planet 
orbiting a moderately bright (V = 11.1), G5 star 
with a super-solar metal abundance ([Fe/H] = 
-1-0.18 ± 0.08). WASP-123 is depleted in lithium 
{log A{Li) < 0.5) suggesting an age of several Gyr. 
This star falls into an area of parameter s pace for 
which gyrochronology is poorly calibrated ( Jeffriesl 
20141) . The Barnes calibration gives an age greater 


than the present age of the universe and Mama- 
jek and Meibom calibrations are not applicable as 
the star’s colour results in the calibrations requir¬ 
ing the logarithm of a negative value. The age 
we derive using BAGEMASS, 6.9 ± 1.4 Gyr, sup¬ 
ports an advanced age. Planets of similar mass 
and radius are not uncommon and are frequently 
found in orbits ~ 3-days around such host stars. 
This makes WASP-123 a typical example of a hot- 
Jupiter system. However, ev en these can prove 
surprising (e.g. WASP-47 : iBecker et al.l 2015t 
Neveu-VanMalle et al. 2016ll and/or contribute as 


vital controls to other studies. 
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Table 5: Stellar and planetary parameters parameters determined from spectra and MCMC analysis. Spectral 
parameters have formal uncertainties while parameters found via MCMC are the median values of the 
posterior distributions with an uncertainty corresponding to the Icr confidence interval. 


Spectroscopic Parameter 

WASP-120 

WASP-122 

WASP-123 

Tycho-2 ID 

8068-01208-1 

7638-00981-1 

7427-00581-1 

USNO-B ID 

0441-0033568 

0475-0113097 

0571-1147509 

RA (J2000) 

04:10:27.85 

07:13:12.34 

19:17:55.04 

Dec (J2000) 

-45:53:53.5 

-42:24:35.1 

-32:51:35.8 

V Magnitude 

11.0 

11.0 

11.1 

Tycho (B—V) colour 

0.523 ± 0.083 

0.78 ±0.11 

0.48 ±0.17 

Spectral Type 

F5 

G4 

G5 

Distance (pc) 

437 ± 21 

266 ± 10 

214 ± 11 

BAGEMASS Age (Gyr) 

2.6 ±0.5 

5.11 ±0.80 

6.9 ± 1.4 

Stellar Effective Temperature, Teff (K) 

6450 ± 120 

5720 ± 130 

5740 ± 130 

Stellar Surface Gravity, loggs 

4.3 ±0.1 

4.3 ±0.1 

4.3 ±0.1 

Stellar Metallicity, [Fe/H] 

-0.05 ±0.07 

0.32 ±0.09 

0.18 ±0.08 

Projected Rot. Vel., V sini (km s”*^) 

15.1 ± 1.2 

3.3 ±0.8 

1.0 ±0.7 

Stellar Lithium Abundance, log A(Li) 

< 1.2 

< 1.0 

< 0.5 

Micro turbulence (km s“^) 

1.5 ±0.1 

0.9 ±0.1 

1.0 ± 0.1 

Macro turbulence (km s“^) 

6.0 ± 0.8 

3.4 ±0.5 

3.4 ±0.5 

MCMG Parameter 

WASP-120 

WASP-122 

WASP-123 

Period, P (d) 

3.6112706 ±0.0000043 

1./ -l-UUOOO_o. 0000026 

2.9776412 ± 0.0000023 

Transit Epoch, To 

6779.43556 ± 0.00051 

6665.22401 ± 0.00021 

6845.17082 ±0.00039 

Transit Duration, ri 4 (d) 

0.1483 ±0.0016 

0.09117 ±0.00082 

0.1289 ±0.0014 

Scaled Semi-major Axis, a/Rs 

5.90 ±0.33 

4.248 ± 0.072 

7.13 ±0.25 

Transit Depth, (Rp/Rs)^ 

0.00655 ± 0.00016 

0.01386 ± 0.00029 

0.01110 ±0.00027 

Impact Parameter, b 

0.78 ±0.02 

0.8622 ± 0.0071 

0.530 ± 0.049 

Orbital Inclination, i (°) 

82.54 ± 0.78 

78.3 ± 0.3 

85.74 ±0.55 

Eccentricity, e 

D np;7+0-022 

U.UO/_o.oi8 

0 (adopted;< 0.08 at 2 ct) 

0 (adopted;< 0.12 at 2 cr) 

Argument of Periastron, cj (°) 

97+48 
^ ' -28 

- 

- 

Systemic Velocity, 7 (kms“^) 

19.836 ±0.013 

34.5934 ± 0.0017 

16.9344 ± 0.0017 

Semi-amplitude, Ki (ms“^) 

509 ± 17 

185.1 ± 2.3 

114.2 ± 2.2 

Semi-major Axis, a (AU) 

0.0514 ±0.0007 

0.03005 ± 0.00031 

0.04263 ± 0.00074 

Stellar Mass, Ms (Mq) 

1.393 ±0.057 

1.239 ±0.039 

1.166 ±0.061 

Stellar Radius, Rs (Rq) 

1.87 ±0.11 

1.52 ±0.03 

1.285 ±0.051 

Stellar Density, ps (p©) 

n 91 9+0.041 

U.Z1Z_Q Q3^ 

0.351 ±0.018 

0.548 ± 0.059 

Stellar Surface Gravity, log((jis) (cgs) 

4.035 ± 0.049 

4.166 ±0.016 

4.286 ± 0.032 

Limb-Darkening Temperature, IpD 

6440 ±120 

5750 ± 120 

5740 ± 130 

Planet Mass, Mp (Mj^p) 

4.85 ±0.21 

1.284 ±0.032 

0.899 ± 0.036 

Planet Radius, Rp (Rjup) 

1.473 ± 0.096 

1.743 ±0.047 

1.318 ±0.065 

Planet Density, pp (pj^p) 

1 r 1+0.33 

- 1 -. 01 _Q 26 

0.243 ±0.019 

0.393 ± 0.056 

Planet Surface Gravity, log(pp) (cgs) 

3.707 ± 0.056 

2.985 ± 0.022 

3.07 ±0.04 

Planet Equilibrium Temperature, Tcq (K) 

1880 ± 70 

1970 ± 50 

1520 ± 50 
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Fig. 10.— Masses, eccentricities and periods 
of transiting planets with non-zero eccentricities 
quoted in literature and masses > O.SMjup mass. 
The most convincing eccentricities are those as¬ 
sociated with more massive planets. Notable ex¬ 
amples a nd WASP-120 b ar e labelled. Data from 
TEPCat ( Southworth 2011 1. 
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